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computer operations for calculating the correlation inte-. 
gral, as compared, with direct methods. . 
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RAPT.D CORRELATION ANALYSIS IN A GROUP OF MOTIONS AND 
HOMOGENEOUS SCALE TRANSFORMS OF THE PLANE (M(2)xR + ) 

D. K. Tkhabisimov 


A method is proposed for picking out a given image in a plane, /2* 
irrespective of translations, rotations, and homogeneous scale 
transforms of the plane. The correlation integral is calculated, 
using the methods of abstract harmonic analysis. In this case, 
and also for groups of translations, it is possible to carry out 
the calculations using the rapid Fourier transform (BPF) . The 
latter fact makes it possible to substantially reduce the number 
of computer operations for calculating the correlation integral, 
as compared with direct methods. 

Introduction /3 

Schwartz's inequality (1] lies at the base of correlation 
analysis. For any functions f(x) and $(x) of L? 

i , . -'O if ex) doc d d-~ J ¥*(*) dx)!^ 

We will assume that the functions f{x) and 4>(x) are determined 
in a homogeneous Euclidean two-dimensional space X, with a group 
of motions G (G is locally-compact) . Then, for any^EEG, Schwartz's 
inequality is also fulfilled: 

J - H: ° ^ ( f‘ x) i%)d* fy 2 (fU) dx ) a) 

X / f 

If the measure dx is invariant relative to the group G, then the 
right-hand portion of (1) does not depend on Insofar as, in 
practical problems, the detection of an invariant measure is 
difficult, we will make use of the Euclidean measure dxidx 2 , 


*Numbers in the margin indicate pagination in the foreign text. 
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relative to the invariant measure in X (2). The problem or' the 
identification of images amounts to a search for that element of 
for which one of the local maxima of the normed correlation 
function is achieved. 

A iy) - 7 , w 

is the Jacobian transforms x^rx. 

The methods of harmonic analysis in groups 12,3] ma’e it 
possible to reduce the calculation of the correlation function 
to integral transforms, in the capacity of the kernel of which 
are matrix elements of the nonreducible unitary concepts of the 
corresponding groups [4]. The number of calculations is reduced 
substantially with numerical solution of problems of identifi- 
cation of the images, based on harmonic analysis. For example, 
for calculation of the correlation function in a group of motions 
of the plane M(2) , with identification of images irrespective of 
their translations and rotations in the plane, N 1 * 3 * log,>N operations /4 
are required, as opposed to^N 5 with calculation of the correlation 
integrals by the known methods [5] (N is the number of elements 
of the discrete network according to each of the parameters) . 

Proposed in the present study is a method of rapid calculation of 
the correlation function in a group of motions and homogeneous 
scale transforms of the Euclidean plane (M(2)xR + ) for identifi- 
cation of the images, irrespective of their translations, rota- 
tions, and homogeneous scale transforms in the plane. Also ex- 
amined are translations and scale transforms of the images in the 
plane . 


where J(0f) 


1. Harmonic Expansion of Correlation Function 

Let X be a homogeneous space with a local ly-compact group 

of motions G. It is common knowledge that X may be realized as 

a space of left classes of contiguity, according to the stationary 
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subgroup H of some point aex [4], In this case, the transforms 
are given by the formula 9* . The 

harmonic expansion of the functions in a homogeneous space X is 
reduced to the expansion of the functions in the group G, which 
are constant in the left classes of contiguity according to the 
subgroup H, i.e., those in which f Jp*-****- 

The correlation function in the group G has the form 


Ktf ) - / muff) A, 


K*>) 




where are the functions which are constant 

in the left classes of contiguity according to the subgroup H, 
and the measure dh is invariant in the group G. We will conduct 
spectral analysis of function (3) . For this purpose, we will 
avail ourselves of the harmonic expansion of the functions f(^) 

, which has the form: 

d* j iu 


and 


//«)•= Z Z *£ 

<k 

- z z 




=i 


/ 

7* 
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where 


(tu - J ♦ 

^ ^ {fi 


(6) 


(?) 


ORIGINAL PAGE IS 
OP POOR QUALITY 



(the line above designates complex conjugation) 


Here is the set of mutually nonequivalent nonreducible 
unitary representations of class I j4J, and 1 Up) are the matrix 
elements of the columns which corresponds to the basis vectors 
so 




that ; d (j is the dimensionality 


of the matrix ||tji f>\\ 


Utilizing (4) - (7), we will calculate the coefficients of 
harmonic expansion of function (3). We will substitute expression 
(3) into the equation 



hm 
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( 6 ) 


in place of the function K(h). By changing the order of integrating 
in (8) and utilizing the expansion 

4 * - £ 4 *, fa). 

we obtain: 



It follows from the unitariness of representation T* > , as well 
as from (6) and (7), that 



Then, the expansion of the correlation function is written as: 
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(XI) 


KQhl ti, £(})■ 


Formula (11) may be written in matrix designations 

K($) = Jz (c* T**), 

where c^KVl, In virtue of the fact that 

the functions f(^D and + (^) and the matrix elements t? \ are 
constant in the left classes of contiguity, according to the 
stationary subgroup H, the integrals (6) and (7) may be rewritten 
in the form of integrals according to the homogeneous space X=G/H 
[4]. Expansion (11) serves as the basis for constructing algo- 
rithms for calculating the correlation function. 

2. Scale Transforms and Translations of Images in the Plane 

Let (X is the Euclidean plane) . 

We will determine the effect of the group G on X in, the following 


manner: 




The normed correlation function in the group G has the form: 

Kfaifati, ’x)f(x) Jx. < I2 > 


The group G, acting on the plane, is the direct sum of the two 
groups of linear transforms of the line; therefore, without 
limiting the continuity, one can study the expansion of the 



correlation function in the group acting on the line. The desig- 
nation for this group, as before, will be G. 


a 


Group G is the hybridized product of the additive group of 
real number (R) and the multiplicative group (R + ) [4]; utilized, 
therefore, for harmonic expansion are the transforms of Fourier 


%, m o, 

and Mellin 

M ati ffi) = * ftj) y/ w 

where Ren«0. Thus, let 

K(a, < f(fx)ffx.> dx , U53 


where ^{a^)-x=axf{ / a 70 • 

(15) to the form 


The Fourier transform reduces 


Kfa, i) - Qty) F fay). 


U6) 


where 


0 ( y ) - ^ ' ^**1 = 4 *. 

In the given multiplicative form, k(a,y) is easily calculated, and, 
having done the inverse Fourier transform, one can obtain (15) 
from its expansion. However, the presence of a scale factor in the 
right-hand part of expansion (16) leads to the necessity of inter- 
polation of the function F(ay), with its determination in a dis- 
crete network. This difficulty is eliminated if one applies the 


6 



Mel 1 in transform to (16). 


We will divide the line into orbits, which are homogeneous 
relative to R + . We will determine the correlation function in /8 

each of the orbits in the following manner: 


a / 

Kfavj « / P-fr, -/*», X 

' L K (e,o) , > m ° 


(i?) 


I# 


The spectral expansion of function (17) on the haif lines y>0 and 
y <0 is calculated using transform (14), and has the form: 


S(0,\»=M, k 





( 18 ) 


for y>0 


sou* W A*? ('V} <“ 


for y<0. Now, it is simple to recreate function (15) from ex- 
pansion (18) - (19) , if one uses the inverse transforms of Mellin 
and Fourier, and if one also notes that K(a,oU ‘ F{0) • T he 
lines utilized here, and the inverse Fourier and Mellin transforms, 
are implemented quickly (in the sense of BPF) [5], which leads to 
a substantial decrease in the number of calculations. This is 
associated with the fact that Mellin' s transform, after substi- 
tution into (14) of^=!L t , is reduced to the Fourier transform of 
the function If one designates the function f (^) in a 

finite set of points N, in the interval (-£ ,«, ), then, selecting 
the points of the readings in the following manner =■ (I-X)AT 

( t • = " Xj A/~J i ^ ), one can implement transform (14), 
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using BPP. 


3. Cor relation Analysis in a Group of Motions and Homogeneous /9 

Scale transforms of the Plane (M (2) xR^ ) ~ 

According to (2), the normed correlation function in the 
group of motions and homogeneous scale transforms of the plane 
(M(2)xR^) has the form: 

K' IWifaft) - ffx) Jx , <20) 

where fhc , x and a are the 

vectors with coordinates xi,X 2 and ai,a?, respectively; x a desig- 
nates the vector rotated at an angle a, R is the coefficient of 
scale transform, and a<. and a z are the translation parameters. 


We will first carry out spectral expansion of function (20) 
in the group M(2) , which is done using the Fourier-Bessel trans- 
form: 

K lf kfi - [iw 

and Fourier's discrete transform: 



As shown in study [5], the Fourier form of function (20) has the 
form: 




where 
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Kfie, f,r /a i) = S(e, ?;*,*), 

«i - pceSfi a : =J>i*r. 

After the application of Moilin' ■ transform to (23), the final 
expansion of the correlation function (20) is written in the 
following manner: 



where 







n 

-U,y 



The integral transforms in (25) , which are used to carry out 
expansion (24), are reduced to BPF through discrete realization 
(see paragraph 2) . If the initial images are given in a discrete 
network N x N in size, and if N points are also taken for each of 
the parameters of the group, then it requires «~N‘'log ; >N operations 
for calculation of the correlation function (20) , as opposed to 
^N ! ’, if the integrals in (20) are calculated for each ^ by the 
direct methods. 

Conclusion 

The use of the methods of abstract harmonic analysis for the 
identification of images on a computer makes it possible to con- 
side 3.v.bly reduce the number of calculations. The methods indicated 



in the present study make it possible to rapidly identify an 
image, irrespective of its rotations, translations, and homo- 
geneous scale transforms in the plane, which may prove useful, 
for example, durirn the reading of print texts on a computer. 

The author expresses his gratitude to B. I. Kolosov and 
D. A. Usikov for their demonstrated interest in the study, and 
their valuable scientific consultation. 
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